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Abstract 

We consider hybrid states with heavy quarks in the frame of re- 
cently proposed constituent gluon model. Limiting regimes for model 
Hamiltonian are discussed. It is shown that these regimes match 
smoothly, hence predictions of the model are definite enough. We pre- 
dict lowest ccg hybrids at 4.1 it 0.1 GeV and bbg hybrids at 10.5 ± 0.1 
GeV. Regge-trajectories for hybrids are derived from the formalism. 
Experimental aspects of our predictions are discussed. 
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One of the main challenges to experimentalists in hadron physics is that 
existence of gluonic degrees of freedom had not been confirmed yet by direct 
experiment in the nonperturbative regime. Presence of gluonic fields in the 
QCD Lagrangian make us think that they should reveal themselves as a 
constituent particles in the form of glueballs and hybrids. The important 
questions are : how can wc distinguish an object containing valcnt gluon 
from ordinary meson and baryon resonances and where can we find this 
object ? Hybrids have more attractive features (with respect to glueballs - 
see [1] ). Most of models predict hybrids with light quarks at 1.5-2 GeV , 
but there are too many resonances at this region and experimental situation 
is far from obvious (despite the fact that standard nonets arc overpopulated 
now). On the other hand , charmonium and bottomonium spectroscopy is 
well understood now and one can identify heavy exotic much easier than the 
light one. These are the reasons for intensive studies of heavy quark hybrids 
carrying now in the literature [2]. 

Recently constituent gluon model based on Vacuum Background Correla- 
tors method [3] and QCD string model [4,5] was proposed [6,7]. This model 
was applied to heavy hybrids spectroscopy in [8] where masses of lowest and 
states were estimated. In this letter we discuss potential and string regimes 
for QCD string in the case of heavy hybrids. 

Consider adiabatic approximation for heavy hybrids. For beginning we 
neglect the Coulomb interaction. Then it is natural that "fast" gluonic sub- 
system energy is of order ^Jo [a is string tension) and quark subsystem energy 

should be ^ times square root from the mass ratio i.e. E^q ~ \f^^ jva 
where m is heavy quark mass. For short distances interquark potential must 
be oscillator one. It's easy to see that unharmonic corrections will be of 
the next order on adiabatic expansion parameter [8]. ( For interquark dis- 
tance p much larger than interquark potential should be crp, of course.) 
Therefore gluon subsystem energy can be obtained after substraction of the 
oscillator interquark potential from the full Hamiltonian. 

Coulomb interaction almost doesn't change nothing. Interquark Coulomb 
repulsion is negligible because its parameter a^/G is very small. Coulomb 
quark-gluon attraction change gluon energy more drastically but still gluon 
subsystem energy is proportional to ^ . 

Hamiltonian of the gluonic subsystem depending on p as a parameter 
(neglecting spin degrees of freedom) has been obtained in [8]. In short, main 
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assumptions made are :1) Quenched approximation for quarks and gluon ; 2) 
Area law for Wilson loops ; 3) Forbidden time backtracking for all particles 
[4]; 4) Straight-line ansatz for minimal surface in the area law [4]; 5) S-wave 
for quarks [8]. The resulting Hamiltonian reads (in Minkowski space-time ) : 

rr _ Pr fig + ai + a2 | ^ , 



12(/x, + 2n)V 6//2r2 



where 



a= f\(5y{(5)-d^ t ^v^= t d(5v{(5){(5 - if 
Jo Jo iy[p) Jo 

Here p is interquark radius ; r is the gluon 3-coordinate with respect to 
quarks c.o.m. (quarks have the same mass) ; fi and f2 are gluon-quark and 
gluon- antiquark 3-coordinates ; (3 parametrize the space part of the minimal 
surface (see [4]) ; z/i,z/2 and /ig are auxiliary fields ( effective string energies 
and gluon mass) ; Lgi is the angular momentum of gluon - it is integral of 
motion in the limit p ^ 0. Strictly say, Lgi is not integral of motion cause 
condition p << r is not satisfied for c- and b-quarks ( they are not enough 
heavy - for our purposes ). We wish to emphasize that the gross structure 
of hybrid spectrum is obtained in the limit p ^ ; p^/r^ terms are cor- 
rections of order 400 MeV for charmonium hybrids and about 200 MeV for 
bottomonium ones. For real hybrids angular momentum is distributed any- 
way between quark and gluon subsystem , but general structure of spectrum 
remains unchanged. 

In the limit p — > one can let ui — V2 and obtain gluon subsystem 
Hamiltonian in zeroth approximation : 

This is a Hamiltonian of heavy-light string system described in [5] with 
effective string tension 2cr - two quark-gluon strings merge into the only 
string. There are two opposite hmits for Hamiltonian (2) described in [4,5]. 
Namely, for Lgi — one can integrate over auxiliary fields [ig and v in the 
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path integral representation for qqg Green function (see [4,5]) and obtain : 

Hgi^yffi + '^^r (3) 

This is a "naive" Hamihonian for zero- mass particle confined by scalar po- 
tential ar (potential limit appeares also for constituents much heavier than 
^/a ). The opposite limit is string limit when r is not dynamical variable - 
this is a limit L » 1. Then one is to replace r by its extremal value : 
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r'=i-\hf. (4) 



and neglect as it is usual [4]. After that one can integrate over field v 
replacing it by its extremal value : 



(5) 



with the result : 



E^pTTa^L{L + l) (6) 

So, as it was shown in [5], string regime for Hamiltonian (2) leads to 
Regge trajectories with the correct slope. For our case with effective string 
tension 2a the slope of Regge trajectories is 2na. Note that result (6) is 
obtained by means of adiabatic approximation, hence it is valid for 

1 << Lgi « m/y/^ (7) 

Condition (7) is satisfied only for bbg hybrids and results (4)- (6) are applicable 
only in this case. We do not consider pure string case with Lgi » 1 when 
even heavy quark mass is negligible with respect to string mass (6). 

These marginal regimes arc purely academical . Indeed, gluon can't have 
Lgi = cause physical states for gluon are electric and magnetic (see [7]). 
On the other hand, here we are interested only with lowest states with < 
L'^i >— 2 and do not evaluate masses of highly excited states. 

In [2] variational procedure for the auxiliary field p was used. In this letter 
we'll show that both potential and string regimes gives the same (within 
accuracy of about 50 MeV) predictions for lowest states. 
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Indeed , variational calculations for Hamiltonian (2) with v = ar gives 
for gluon subsystem energy (without Coulomb forces) 

Egi = 4.36v^ fa 1.7GeV (8) 

String energy (6) for Lgi = 1 is 

E^f = v^V27r\/2 = 2.98^/CT 1.2GeV (9) 

is relatively small with respect to (8). Actually correction to regime (6) is not 
small. Following the method described in [4] one can estimate first correction 
to the string regime: 

SEgi = v^^l^V^ ^ l.SSV^ (10) 

Adding (9) and (10) one obtains gluon subsystem energy calculated by means 
of string regime formulaes: 

Egi = 4.51v^ fa l.SGeV (11) 

As it is evident from (8) and (11) potential and string regimes match smoothly 
and gives the same predictions for gluon energy of lowest hybrids. This is 
a rather general property of the QCD string formalism developed in [4,5] . 
It's instructive to estimate the deviations from regime (6) because of quark 
correction. According to (1) interquark potential arise from the terms: 

, nil np'\ 
'''^'^''US2(,, + 2n)V-;SF/., '''' 

9'' 

after averaging over gluon wave function obtained from Hamiltonian (2). It's 
easy to see that for regime (6) : 

d-l/^/L; (13) 
For the last term in the straight-line ansatz [4] one can estimate : 

- 1/^ (14) 
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Therefore the last term in (12) decreases most rapidly and there is no effective 
interquark repulsion (cause there is no repulsion for Lgi = 1). Also it's evident 
that quark correction to regime (7) decreases as l/v^ (i.e. very slowly) and 
for regime (7) quark correction to the gluon energy is about 100 MeV. 
Three Coulomb terms should be added to the full Hamiltonian: 

Vc„, = (i + i) + (15) 

2 Vri r2/ 6 p 

So,to take into account Coulomb forces one is to add the term —Sag/r to 
zeroth-approximation Hamiltonian (2). The explicit form of quark subsystem 
Hamiltonian is: 

H,,-^ + Ap' + ^^ (16) 
m op 

where oscillator term of the potential is derived in (12). We calculate lowest 
eigenvalues of Hamiltonian (16) variationally and designate them as Egg . For 
the calculation of the absolute value of hybrid mass we use the procedure 
proposed in [7,8] : the additive constant in the potential is twice larger 
than for heavy-light system. Let us define the constant term from masses 
of P-states of heavy-light system . The Hamiltonian for heavy-light system 
coincides with hybrid Hamiltonian (2) after replacing 2a — > a. Let us neglect 
Coulomb for beginning. Since energy Egi is proportional to ^/a, heavy-light 
gluon energy is y/2 smaller than for hybrid system . Therefore: 

Mggg ^ 2Mgg{L = l) + Egl{l " ^) + Egg (l7) 

Here Egi is gluon subsystem energy discussed above , Mgg is the experimen- 
tal value of heavy-light state mass with L = 1 (2.42 GeV for D-meson ) . 
Actually heavy-light and hybrid Hamiltonians differ by Coulomb terms and 
one is to use a function Egi{as) ■ Then 

Mggg = 2Mgg + («,) - ^Egl (4a,/9) + Egg ( 18) 

Prom (18) it's clear that predicted hybrid mass in our model depends on 
quark mass very slowly cause the only mass dependence is contained only in 
the last term in (18) and this term is inversly proportional to -^/m . Numerical 
results for charmonium hybrids are listed in the Table I. 
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Table I 

Hybrid masses calculated by means of string and potential regimes 
formulaes. The parameters are : cr = 0.18 GeV^, as — 0.3 

Quark mass,GeV Potential, GeV String,GeV 
rUc = 1.2 4.209 4.134 

1.3 4.198 4.125 

1.5 4.178 4.117 

1.7 4.162 4.100 

So, for charmonium hybrid mass we anticipate : 

Mccg = 4.1 ± 0.1 GeV (19) 

Since P-levels for B- mesons are still unknown wc can't use (18) for evaluation 
of bottomonium hybrid mass. Therefore we are to use the same procedure 
as in [8] i.e. define additive constant from S-levels of B-mesons. Hence, our 
predictions for ccg mass coincide with predictions of [8] . 

We wish to discuss some aspects of flux-tube model [9] in comparison 
with our results. First, there is no Regge-trajectories in flux-tube model be- 
cause it is purely nonrelativistic . Second , in the original papers [9] small 
oscillations approximation for string excitations was used. Numerical cal- 
culations in the frame of "one-bead flux-tube model" [10] show that this is 
not the case. Moreover, rather general arguments given in [8] prove that for 
(infinitely) heavy quarks gluon subsystem oscillation is much larger than for 
quark subsystem. It's sufficient that small string oscillations never appear in 
the flux-tube model ; hence , the approximation made in [9] is rather doubt- 
ful. Our predictions for heavy hybrid mass coincide with flux-tube model [9] 
( and with lattice calculations [11]) , but it follows from the fact that various 
prescriptions for additive constant were used . 

Quantum numbers of our hybrids are (see [7]): 

jpc_ QT+^1T+^ 2^+, 1^- (20) 

Our hybrids can be produced directly in e"^e~ collisions ( J^^ — 1 is al- 
lowed ) , but appropriate width should be very small cause hybrid itself could 
be produced through mixing between ordinary and hybrid meson states . 
Of course, "gluish" channels like pp annihilation and highly excited quarko- 
nium decays are preferrable. Hybrid states should be very narrow cause for 
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ccg hybrids their masses are close to open charm thresholds and decays are 
strongly supressed by phase space ; for bbg hybrids decays to open bottom 
are forbidden (strongly suppressed , at least). 
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